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Abstract 

Recent results from T2K, MINOS and Double CHOOZ all indicate a sizeable reactor 
angle ^13 which would rule out conventional tri-bimaximal lepton mixing. However, 
it is possible to maintain the tri-bimaximal solar and atmospheric mixing angle 
predictions, 9i2 ~ 35°, ^23 ~ 45° even for a quite sizeable reactor angle 6*13 8°, 
using an ansatz called tri-bimaximal-reactor (TBR) mixing proposed by one of us 
some time ago. We propose an explicit ^4 model of leptons based on the type I 
seesaw mechanism at both the effective and the renormalisable level which, together 
with vacuum alignment, leads to surprisingly accurate TBR neutrino mixing, with 
the second order corrections to mixing angles having small coefficients. 
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1 Introduction 



Recently T2K have published evidence for a large non-zero reactor angle [T] which, when 
combined with data from MINOS and other experiments in a global fit yields [2] , 

013 = 8° ±1.5°, (1.1) 

where the reactor angle is defined in the usual PDG convention [3] and the errors indi- 
cate the one a range, although the statistical significance of a non-zero reactor angle is 
about 3(7. A non-zero reactor angle is also consistent with the first results from Double 
CHOOZ [4J. 

If confirmed, these recent observations would rule out con?;ent«ona/ tri-bimaximal (TB) 
mixing which predicts a zero reactor angle [5]. Many alternative proposals [B] have been 
put forward since the first T2K results, all aiming to accommodate a non-zero 6'i3. Typ- 
ically, such models also yield (sizable) deviations to the other predictions of TB mixing, 
namely tri- maximal solar mixing 612 = 35.26° and maximal atmospheric mixing 623 = 45°. 
However, since the latter remain in good agreement with current global fits [2], it would 
be nice if, somehow, these good predictions of TB mixing could be maintained while at 
the same time allowing a non-zero reactor angle. In fact the idea of maintaining the TB 
predictions for the solar and atmospheric angles, while switching on the reactor angle, 
was suggested some time ago (before T2K results) by one of us and was referred to as 
tri-bimaximal- reactor (TBR) mixing [7]. 

The TBR ansatz [7] postulates a free reactor angle 613 but with fixed S12 = l/y/S 
and S23 = corresponding to tri-maximal solar mixing 612 = 35.26° and maximal 

atmospheric mixing 623 = 45°. In the PDG convention for the PMNS mixing matrix |3] 
the TBR ansatz for the mixing matrix is then [7j, 



TBR 



-^(1 + 72.136^^) ^(1-^^.136^^) >3 

V ^(i-V2^i3e^^) -73(1 + 75^136^') >3 y 



;i.2) 



■ "^21 ■ '-^31 \ m 

where P = diag(l, e*~, e*~) contains the usual Majorana phases. Note that TBR 
mixing reduces to TB mixing in the limit that 613 — )■ 0. However in general TBR mixing 
involves an arbitrary reactor angle 613 which could in principle be large, without causing 
any deviations from TB solar and atmospheric mixing S12 = l/\/3 and S23 = 1/v^- 

The TBR ansatz in Eq. (11. 2p is clearly very simple to write down. The obvious 
question is whether there is any model that can give rise to TBR mixing. Such a model 
of TBR mixing should also explain the smallness of the reactor angle ^13 as compared 
to the tri-maximal solar mixing angle 612 = 35.26° and the maximal atmospheric mixing 
angle 6'23 = 45°. At the same time that the TBR ansatz was proposed, a mechanism was 
suggested called partially constrained sequential dominance (PCSD) which could lead 
to TBR mixing [7j. Although no actual model was proposed, it was shown how PCSD 
could originate due to a distortion of the vacuum alignment which was previously used to 
account for TB mixing via constrained sequential dominance (CSD), which in turn could 
originate from family symmetry [7]. 
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In this paper we propose the first exphcit model of leptons based on the type I 
seesaw mechanism at both the effective and the renormahsable level which, together 
with vacuum alignment, leads to PCSD and hence TBR mixing. To understand the 
approach we are following in this paper, it is useful to first recall the CSD approach to 
TB mixing. The basic ingredients there are sequential dominance (SD) [8l|9] and vacuum 
alignment. The strategy of combining SD with vacuum alignment is familiar from the 
CSD approach to TB mixing [TU] where a neutrino mass hierarchy is assumed and the 
dominant and sub dominant fiavons responsible for the atmospheric and solar neutrino 
masses are aligned in the directions of the second and third columns of the TB mixing 
matrix, namely (1, 1, —1)^ and (0, 1, 1)"^. The idea of PCSD [7J is to simply maintain the 
subdominant fiavon alignment, ipp^ = (1, 1, —1)"^, while considering a small perturbation 
e to the dominant fiavon alignment, (p^,^ = {e, 1, 1)"^. Assuming PCSD, with a strong 
neutrino mass hierarchy, Im^l > |m2| ^ ^ 0, leads to a perturbed neutrino mass 
matrix rrii, [7], 

= ^(f),y2(j)^J' + ^(pus^'uj', (1.3) 

where and are the leading order neutrino mass eigenvalues. The explicit ^4 models 
we propose involve the accurate fiavon alignments (p^^ = (1, 1, —1)"^ and (j)^^ = {e, 1, 1)"^ 
which appear quadratically and reproduce the neutrino mass matrix as in Eq. fll.3p quite 
accurately in the case of the renormahsable model. 

Another important result of our paper concerns the analytic diagonalisation of the 
neutrino mass matrix in Eq. (11. Sp to second order e^. Previously it has been shown that, 
to leading order in |m2|/|m3| and e, Eq. (ll.3p leads to TBR mixing in Eq. (II. 2p [7|[TT]. 
Remarkably, we find that the coefficients of the second order corrections to the mixing 
angles are suppressed, making the approximate TBR mixing resulting from Eq. fll.3p to 
be more accurate than expected. Thus the renormahsable A4 model predicts surprisingly 
accurate TBR mixing with 612 ~ 35°, 623 ~ 45° even for quite sizeable ^13. 

We note that PCSD is based on the general approach to model building known as 
the indirect approach [12] which involves the quadratic appearance of the fiavons as in 
Eq. fll.3l) . Recently an alternative type of model of TBR mixing has been proposed based 
on 5*4 family symmetry, with a type I plus type II seesaw mechanism, although no detailed 
model was proposed and the necessary vacuum alignment was not studied [13] . Moreover, 
TBR mixing does not follow in a general way from the model in [13], but only occurs for 
special choices of parameters. 

The layout of the rest of the paper is as follows. In Section [2] we propose an explicit 
A4 model of leptons at both the effective and the renormahsable level which, together 
with vacuum alignment, leads to PCSD and hence TBR mixing. In Section [3] we give the 
results of the analytic diagonalisation of the neutrino mass matrix in Eq. (I1.3p to second 
order e"^. The straightforward generalisation to the case with a non-zero lightest neutrino 
mass mi 7^ is presented in Appendix |X] and might be relevant for some future model. 
Section H] concludes the paper. 
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Table 1: Lepton, Higgs and flavon fields in an indirect model. 

2 A4 models of tri-bimaximal-reactor mixing 

In this section we present the details of the model, first at the effective, then at the 
renormalisable level. As pointed out in the introduction, our models follow the indirect 
approach and the aligned flavons appear quadratically in the light neutrino mass ma- 
trix vHy as in Eq.f ll.3p . In order to achieve this via the type I see-saw mechanism, and 
to ensure a diagonal charged lepton mass matrix, additional Z4 shaping symmetries are 
also employed. In such models, the discrete family symmetry is pivotal in two aspects: it 
ii) helps generate the required flavon alignments and iii) combines the three generations 
into multiplets of the group so that, together with the corresponding product rules, one 
obtains mass matrices which depend on only a few free parameters. In this paper we 
choose to work with ^4 as it is the smallest non-Abelian finite group with an irreducible 
triplet representation. We apply the ^4 basis in which the triplets are explicitly real as 
given for example in [H]. Denoting a general A4 triplet as c = (01,02,03)^ and defining 
ijj = e^'^'/^, the product rules can be summarised as 

+ cA , (2.1) 

\C2C'l/ 

corresponding to two triplets and the sum of the three one-dimensional irreducible repre- 
sentations Ir, with lo = 1 being the trivial singlet. Furthermore, Ir Is = l(r+s)mod3- 

2.1 The effective A4 model 

The complete list of lepton, Higgs and flavon fields introduced in our model is given in 
Table [H Similar to the A4 model of [15] we have a ?7(1)_r symmetry as well as several Z^ 
shaping symmetries. 



El 

r=0 



C^C^ + U ''C2C2 + W^CgCg) + 
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The neutrino part of the effective superpotential reads, 

W-f ~ Y.iLH^^Nt + mm'^), (2.2) 

where the effectively allowed mixing term N2N^Lpy^ipy^ must be forbidden by the choice of 
appropriate messengers. Notice that there are only two right-handed neutrinos living in 
the singlet representation of ^4. Hence, the model will feature one massless light neutrino. 
Inserting the flavon VEVs 




whose alignment is discussed later, leads to the following Dirac and right-handed Majorana 
neutrino mass matrices 

M3 j ' ^ > 





f 0-2 


0-3 




rriD = 




0-3 


1 , Mr 




{-a2 


"3 ) 





where Oj ~ and Mj ~ xp-- Note that Mr is diagonal by construction. 

Using the type I seesaw formula we can express the light neutrino mass matrix as 

,„„ = ™,M-™5 = ^ _l j + _i 1^ J 1 1 j . (2.5) 

With this structure we arrive at TBR mixing in the neutrino sector as will be shown 
analytically in Section [31 

In order to account for the charged lepton mass hierarchy we identify the right-handed 
charged leptons with ^4 singlets, and distinguish them using the Z4 shaping symmetries. 
The resulting effective charged lepton superpotential then takes the form 

<f ~ -^H,{L^rT' + L^,^^ + L^,e'') . (2.6) 

Inserting the fiavon VEVs 

(<^r) = Vr , {V,) = V, , (<^e) = V, ^ , (2.7) 

whose alignment is discussed later leads to 

<f ~ -^H,{v,L^T' + v^L2^l^ + v,L^e') , (2.8) 

thus yielding a diagonal charged lepton mass matrix m^. In this model, the hierarchy in 
the charged leptons remains unaccounted for, however it is straightforward to implement 
the Froggatt-Nielsen mechanism [TB] to cure this. For the purpose of clarity we will ignore 
this issue in the following. 
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Table 2: Messenger fields of the renormalisable A4 model. 



2.2 The renormalisable A4 model 

As emphasised in pTTj, any non-renormalisable operator of an effective superpotential 
should be understood in terms of a more fundamental underlying renormalisable theory. 
Without such a UV completion of a model, higher order terms which are allowed by the 
symmetries may or may not be present. Thus a purely effective formulation would leave 
room for different physical predictions. We have constructed a fully renormalisable theory 
of the lepton sector. The required messengers are listed in Tabled while the driving fields 
which control the alignment of the flavons are presented in Table [31 

With the particle content and the symmetries specified in Tables [T]l2l we can replace 
the effective neutrino superpotential in Eq. (12. 2 p by a renormalisable one which includes 
the messenger fields, 

3 

i=2 

3 

+ Y.^M^^^^x'i + Mrj:T^) . (2.9) 

1=2 

Integrating out the messenger pairs Xj, Xi and T^, T^, the effective operators of Eq. (12. 2p 
are uniquely generated. Notice that the messenger pairs Tj,T^ do not lead to the afore- 
mentioned mixing term A^IA^Iy^i^j^jj^g. 

The charged lepton sector is formulated at the renormalisable level using only one new 
pair of messengers, fl and Q'^. With the particles and symmetries listed in Tables [1] and [2] 
we get the renormalisable superpotential for the charged leptons 

W^A^ ~ LHan + nVrr' + l^V^/u' + f^VeC' + Mn^n" , (2.10) 

where we have suppressed all order one coupling constants. Integrating out the messen- 
gers, we are led uniquely to I^^'f of Eq. fl^ . 
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Table 3: Driving fields of the ^4 vacuum alignment. 



2.3 Vacuum alignment 

So far we have only postulated the particular alignments of the neutrino-type flavons, given 
in Eq. fl2.3p . and the flavons of the charged lepton sector, Eq. (12. 7p . In this subsection we 
explore the driving potential and prove that the assumed fiavon alignments can in fact be 
obtained in a relatively simple and elegant way. 

The renormalisable superpotential involving the driving fields necessary for aligning 
the neutrino-type flavons is given as 

^flavon.i. ^ A^^{gi'^^^^^^ + g2^y2i) (2.11) 
+ Oau^gr'^dVus 

+ D{g8(Pe^U3 + g9Vr^) + ^'g^VvaV^^^a + g^^X'^'uz + MsSS^ . 

The first line of Eq. (12. lip produces the vacuum alignment (v^j^a) oc (1, 1, —1)"^ of Eq. (12.31) 
as can be seen from the F-term condition^ 

({Vu2)2{'{>U2)z\ f{'^U2)l\ M 

2^71 (^..)3(<^..)i +g2{0 {^^2)2 = . (2.12) 
The terms in the second line of Eq. (12.111) give rise to orthogonality conditions which 



^We remark that the general alignment derived from these _F-term conditions is {(pi^2) (^1; ^1)"^- 
One can, however, show that all of them are equivalent as the resulting mass matrices will be related 
by appropriate changes of the unphysical phases in the matrix P' . Our choice leads to 6e.^i,T all being 
zero in the limit of vanishing e, see Eqs. p. 13113. 15)) . 
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uniquely fix the alignments of the auxiliary flavon fields if a and <^a; 

i^,)^ ■ (ipa) = {^u,V ■ {Va) = ^ {^a) « ( 1 I , (2.13) 




Here we have assumed the fiavons and (f^,^ to be already aligned as in Eqs. (12. 7p and 
(12. 3p . respectively. Finally, the neutrino-type flavon ipi,^ gets aligned by the remaining 
terms of Eq. (12. lip . A vanishing F-term of the driving field Oauj requires 

(^.a) = , (2.15) 

where ni and n2 are independent parameters. They are further constrained by the F-term 
condition of the driving field D which - after integrating out the messenger pair E, S"^ 
and inserting the flavon VEVs - reads, 

gsVeni- ^g9Vr{ggV^2+ 99{0h2 = 0. (2.16) 

Ms 

This shows that rii is naturally suppressed compared to n2 due to the messenger mass. 
With e = ^, we get 

i^u,) = V,, l^lj , (2.17) 

as anticipated in Eq. (12. 3p . 

Turning to the flavon alignment of the charged lepton sector, the renormalisable driv- 
ing superpotential takes the form 

Here we have suppressed all order one coefficients as they are completely irrelevant. 

The triplet driving flelds ^e,/^,r give rise to flavon alignments (v^e.^i.r) with two zero 
components. The singlet driving flelds in turn require orthogonality among the three 
flavon VEVs so that we end up with the vacuum structure of Eq. (12. 7p , possibly requiring 
some redeflnition of the family indices. 

3 Analytic diagonalisation of the PCSD neutrino mass 
matrix to second order 

In this section we diagonalise the effective neutrino mass matrix of Eq. (II. 3p . 

= ^ I 1 1 -1 I I e 1 1 I , (3.1) 
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which arises in the indirect A4 models of Section [21 see Eq. fl2.5p . Numerically the complex 
parameter 

e = ee'^" , (3.2) 

e, (5° G M, is assumed to take an absolute value of e ~ 0.2. It is therefore reasonable to 
expand the mixing matrix which diagonalises m,^ in powers of e. In the limit of vanishing e, 
the absolute values of the mass parameters 

m° = |m°|e-^"° , (3.3) 

correspond to the eigenvalues of m^, and the mixing matrix is exactly of tri-bimaximal 
form. Switching on e changes both the masses and the mixing. In what follows we assume 
a normal hierarchical neutrino mass spectrum which allows us to parameterise the masses 
in Eq. as 

|m°| = ke\ml\ , (3.4) 

with G M being an order one coefficient. The neutrino mass matrix mi, of Eq. (13. ip is 
now diagonalised by the unitary mixing matrix U = P' f/pMNS such thato 

P' = diag(e*^% e*^f , e*"^^) is an unphysical phase matrix which is required to bring the 
PMNS matrix into PDG form, 

/I \ / ci3 si3e-'^\ / C12 si2 0\ 

?7PMNS = C23 S23 1 ~Su OP, (3.6) 

\0 -S23 C23/ V-sise'^ ci3 / V 1/ 

with Cij = cos 6ij and Sij = sin 6ij. The Majorana phases are included in the matrix 
P = diag(l, e*~, e* 2 ). 

We have worked out the mixing matrix U to second order in e using the above notation. 
The results read 

^12 = arcsin^-^, (3.7) 
O23 = j + ik cos(a° -al + 5°) , (3.8) 

= 72 + 573 ^ - «3 + 25°) , (3.9) 

S = f A;sin(a°-a° + 2(5°) , (3.10) 

"21 = "2 ? (3-11) 

asi = «3, (3.12) 

Se = , (3.13) 

(5^ = f A;sin(a°-a° + 5°) , (3.14) 

6r = -4A;sin(a°-a° + (5°) , (3.15) 



^Here we assume the convention in which the effective hght neutrino mass matrix rrii, is defined through 
the bihnear vli^l coupling, in contrast to the convention Vi/Pl adopted in [7]. 
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Figure 1: The change of the solar, atmospheric and reactor mixing angles with increasing e in PCSD. 
This plot, which first appeared without the phase dependence in [9], was obtained using the Mixing 
Parameter Tools of the REAP Mathematica package [TH]. The numerical results shown are consistent 
with our analytic second order equations as discussed in the text. 

leading to a diagonalised mass matrix of the form 

m^'^^ = diag (mi , m2 , ma) = diag ^0 , e , 1 + ^ j l^^al • (3.16) 

Notice that the Dirac CP phase 6 is only given to first order in e as it always appears 
together with sin ^13 whose leading term is already suppressed by one power of e. The 
Majorana phases on the other hand do not receive any corrections of order e or e^. With 
mi = there is only one physical Majorana phase, 

^23 = "21 - "31 = "2 ~ • (3-17) 

To first order in e the deviations in the solar and atmospheric mixing vanish, and only the 
reactor angle ^13 is switched on. Hence the mass matrix of Eq. (13. ip gives rise to TBR 
mixing at leading order. In more detail, the second order results show that the solar mixing 
angle is corrected to values slightly smaller than the TB value, while the atmospheric 
angle can deviate to larger or smaller values, depending on the phases. Furthermore, the 
deviation of 6*23 is bigger by a factor of about 2\/2 compared to the deviation of 6'i2. 

Our analytic expressions confirm the numerical results for PCSD obtained with the 
Mixing Parameter Tools of the REAP Mathematica package [18]. The resulting variation 
of the mixing angles with e is shown in Fig. [H The broadening of the allowed regions 
of the atmospheric and reactor mixing angles is caused by the phase dependence of 6*23 
and ^13, cf. Eqs. (I3.8f3.9p . In the case of the solar mixing angle, our analytic formula, 
Eq. (13. 7p . does not show any phase dependence to second order in e, so the significantly 
smaller broadening of the allowed region of 6'i2 is an effect of third order in e. 
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We can express our results in terms of the deviation parameters s, a and r as defined 

in 



sin0i2 = ^(1 + s) , sin023 = + a) , sin 6^13 = — . (3.18) 

v3 v2 v2 

Comparing these definitions with our second order expressions in Eqs. (13. 7113. 12]) we find 

2 2 2 

s = — ^ , a = ^ k cos(q;23 + 5), r = e + yA; cos(a23 + 25) , (3.19) 

where we have replaced the original phase parameters by the two physical phases 6 and 
a23, keeping only terms up to second order in e. These relations can be combined to yield 
the second order sum rules 

s = -i , a= ^ fccos(a23 + 5) = § ^ cos(a23 + S) . (3.20) 

Using these sum rules and the second order expansion of the PMNS matrix given in [T9] . 
we can rewrite the mixing matrix of PCSD in terms of the reactor deviation parameter r 
and the parameter k = k cos(a23 + S), 



PMNS 



I 



1-V 



re 



r 

u 



V3 

) (1 + 16^-5 + ii^r^ 



J_ / 1 _ 5rf 



r iS 
AS 



1-4k ^2 
12 
-+4> 
12 



) *(1 



3-4k^2 



12 



12 



r 



P. 



(3.21) 

The comparison with the second order expansion of the TBR mixing matrix in Eq. (11. 2p . 
where s = a = 0, 



U' 



TBR 



v 



3(l-f) 

75(1 -re^') 



J_ /1 _ r2\ j;_ 

^/3 4 ; ^/2 



- — fl + -p*' 



J_p-i5 
V2^ 



\ 



P 



(3.22) 



illustrates how accurately TBR mixing is achieved in models of PCSD. 

It is also worth noting that the columns of the Dirac mass matrix in Eq. (12.41) are 
not quite proportional to the columns of the PMNS matrix in Eq. (l3.2ip . so that form 
dominance [20] is violated at order 0{r). This implies that leptogenesis is non-zero [2T] . 
even in the absence of renormalisation group corrections [22j. 



4 Conclusions 

Recent results from T2K, MINOS and Double CHOOZ all indicate a sizeable reactor 
angle 6'i3 which would rule out conventional tri-bimaximal lepton mixing. However, it is 
possible to maintain the tri-bimaximal solar and atmospheric mixing angle predictions, 
6*12 ~ 35°, 6*23 ~ 45° even for a quite sizeable reactor angle such as 6^13 ^ 8°, using an 
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ansatz called tri-bimaximal-reactor (TBR) mixing in Eq. (11. 2p proposed by one of us some 
time ago, along with the neutrino mass matrix in Eq. fll.3p arising from PCSD [7]. 

In this paper we have proposed the first explicit A4 model of leptons based on the type I 
seesaw mechanism at both the effective and the renormalisable level which, together with 
vacuum alignment, leads to the desired form of neutrino mass matrix. After performing 
an analytic diagonalisation of the neutrino mass matrix to second order, we find that the 
coefficients of the second order terms are suppressed, making TBR mixing surprisingly 
accurate for the renormalisable A4 model. The analytic results are confirmed by the 
numerical results in Fig. [1] which illustrates the stability of the atmospheric and solar 
angles as the reactor angle is switched on. 

It is worth emphasising that the A4 models are indirect models involving the quadratic 
appearance of misaligned flavons as in Eq. fll.3l) . Moreover there is no simple Klein 
symmetry that is respected by this mass matrix; in fact, formally, the generators of the 
Klein symmetry involve the real order one parameter k = ^^ly of Eq. (13. 4p , thus proving 
the absence of a direct link between the family symmetry and the Klein symmetry of 
the neutrino mass matrix. Nevertheless, the alignments of the flavons can readily originate 
from a simple discrete symmetry such as ^4, which is broken in a rather complicated way 
in the neutrino sector, although the charged lepton sector is diagonal in the model of 
leptons considered here. 

In a more complete A4 family unified model, for example based on SU{5), one would 
expect the charged lepton sector to be related to the off-diagonal quark mass matrices, 
resulting in additional charged lepton corrections to the lepton mixing angles, together 
with renormalisation group (RG) corrections, as discussed in [23]. However it is worth 
emphasising that, typically, such corrections to mixing angles are not expected to exceed 
about 3° from charged lepton corrections and about 1° from RG corrections [23]. Such 
corrections are not sufficient to account for the observed reactor angle, although they 
may affect the predictions of TBR mixing discussed here, which strictly apply to the 
unrenormalised neutrino mixing angles only. 
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Appendix 



A Analytic diagonalisation of the PCSD neutrino mass 
matrix with to second order 

In this appendix we present the resuhs of the mixing angles and phases for the more 
general case with three non-zero mass parameters m°. The starting point is the PCSD 
form of neutrino mass matrix with non- vanishing m^, namely, 








e 










l+f 1 




1 


i) 












u 


1 





= -i -2 1 -1 + 1 1 _i + U 1 1 . (A.l) 



Clearly, in the limit m\ = 0, this reduces to the PCSD form of neutrino mass matrix 
considered in Eq. (11.31) and Section [31 Adopting the notation of Section [3] and assuming 
the hierarchy 

\m1\ = kie\ml\ , \ml\ = fc2e|"^3l, (A. 2) 

with fcj G M being coefficients of order one or smaller, we obtain to second order in e, 

= arcsm^- e , (A.3) 

= f + |[A;2Cos(a°-a° + 5°)-fcicos(a?-a° + 5°)]e^ (A.4) 

^13 = 73 + 5^ [k2 cos(a° - a° + 25°) + 2h cos(a? - a° + 25°)] , (A.5) 

S = 5°-l[A:2sin(a°-a° + 25°) + 2A;isin(a?-a° + 2(5°)]e, (A.6) 

asi = «3-«i+ 3(fci-4) ^ ' 

X _ gg I fc2fcisin(a°-a;) 2 (' A Q^l 

- "T + 6(fci-fc^) — ^ ' l^-y^ 
5^ = f + sin(a° - + 6') - ^^^^ 1^4-"°^ - h sin(«; - «° + ^°)] 6^ (A.IO) 

5, = f - I [fc2 sin(a° - al + + '^'^g^J'""^ - fci sin(a° - a° + 6')] e\ (A.ll) 

The resulting neutrino masses read 

m^'^^ = diag (mi , m2 , ms) = diag e , A;2 e , 1 + l^^sl ■ (A-12) 

Note that mi can be of same order in e as m2 without changing the mixing matrix to 
first order in e. That is even with ki being of order one, we still find TBR mixing. 
Considering, however, the case where \mW is additionally suppressed by one or more 
powers of e, i.e. ki ~ 0{< e), the above expressions simplify considerably. Then, in 
fact, the ki dependence drops out completely to second order in e, and we recover the 
expressions of Section [3], with /c2 corresponding to k and set to zero. 
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